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Abstract 

We prove that Pade approximants yield increasingly accurate predictions of higher-order coeffi- 
cients in QCD perturbation series whose high-order behaviour is governed by a renormalon. We also 
prove that this convergence is accelerated if the perturbative series is Borel transformed. We apply 
Pade approximants and Borel transforms to the known perturbative coefficients for the Bjorken sum 
rule. The Pade approximants reduce considerably the renormalization-scale dependence of the per- 
turbative correction to the Bjorken sum rule. We argue that the known perturbative series is already 
dominated by an infra-red renormalon, whose residue we extract and compare with QCD sum-rule 
estimates of higher-twist effects. We use the experimental data on the Bjorken sum rule to extract 
a s (M|) = 0.116^QQQg, including theoretical errors due to the finite order of available perturbative 
QCD calculations, renormalization-scale dependence and higher-twist effects. 
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1 Introduction 



Everybody interested in more precise quantitative tests of QCD, or in its place in some Grand 
Unified Theory, would welcome a more precise determination of the strong coupling strength 
a s in some well-defined renormalization prescription, say MS, at some reference energy scale, 
say M Z - Such determinations are normally made using perturbative QCD to interpret data, 
though lattice QCD may also become competitive once systematic effects are better controlled. 
Obtaining the desired level of precision using perturbative QCD requires calculations beyond 
the next-to-leading order. Several processes are calculated to high order in perturbative QCD 
[0,0. However, progress in the high-precision determination of a s [Mz) is hampered by the 
fact that the QCD perturbation series is expected to be asymptotic: 

S{x) = £ c n x n , x = -1 , cn ~ n\K n tf (1) 

n=0 * 

for some coefficients K, 7 |J. Under these circumstances, how can one extract the most infor- 
mation from the QCD perturbation series, and obtain the best value of a s = Ttxl The usual 
answer is to calculate up to order n opt : A n = \c n x n \ is minimized, and use the magnitude A nojrt 
of this minimum term as an estimator for the residual uncertainty. 

^opt 

In this paper we study whether it is possible to estimate S opt (x) = Yl c n x n reliably without 

n=0 

the labour of calculating all the perturbative coefficients c n : n < n opt , and comment on the 
possibility of summing over the higher-order terms S asy (x) = c nX n sufficiently reliably to 
reduce the magnitude of the residual uncertainty below A nopt . Our approach to these issues is 
based on Pade approximants, which are widely appreciated in other areas of physics, and which 
we have recently shown Q can be used to predict higher-order perturbative QCD coefficients 
in agreement with exact calculations (where available) and with the effective charge method 

In this paper, we present new results on the rate of convergence of Pade approximants for 
series of the form ([]]) expected in QCD. We also demonstrate that they reduce significantly 
the renormalization-scale dependence of the perturbative series for the Bjorken sum rule, and 
summarize a comparison with another technique for treating higher-order effects in perturba- 
tive QCD 0. To go further, we transform to the Borel plane, where behaviours of the type 
([!]) correspond to discrete renormalon singularities The Pade technique is a priori well 
adapted to locating such singularities, and we indeed prove that the convergence of the Pade 
approximants is accelerated for the Borel transform of a series such as (|l]). We apply this com- 
bined Pade-Borel technique to the calculated QCD perturbation series for the Bjorken sum 
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rule, and show that it yields a leading infra-red renormalon pole close to the expected location 
in the Borel plane. Assuming this location, we extract its pole residue and use it to evaluate 
the possible magnitude of the infra-red renormalon ambiguity in the perturbative contribu- 
tion to the Bjorken sum rule, which we argue is canceled by a corresponding ambiguity in the 
non-perturbative contribution. We use Pade summation, extracting a s (M|) = 0.116lo;oog from 
the available polarized structure function measurements, including theoretical errors associated 
with renormalization-scale dependence and higher-twist effects. The accuracy of this result 
testifies to the utility of both Pade Approximants and the polarized structure function data. 



Convergence" of Pade Approximants 



We denote Pade Approximants (PA's) to a generic perturbative QCD series S(x) = c n x n by 

n=0 

[N/M] = Y^g^fjg^ : [N/M] = S + O(x^) (2) 

i.e. the PA's are constructed so that their Taylor expansion up to and including order N+M 
is identical to the original series. We have previously pointed out that the next term in the 
Taylor expansion of a [N/M] PA typically provides increasingly accurate estimate c e ^ +M+1 of 
the next higher-order perturbative coefficient cn+m+i of the original series. In the following 
we refer to such estimates as Pade Approximant Predictions (PAP's). 

Let us briefly restate the condition |§ for the convergence of PAP's. With f(n) = lnc n and 
g(n) = d 2 f(n)/dn 2 (where the derivative with respect to n is to be understood in a discrete 
sense), a sufficient condition for PAP convergence, c^ +M+1 — > c^+m+Ii is nm 9( n ) — 0. To 
quantify the rate of convergence, we introduce the quantity 

Cn = C J1^±1 _ 1 = e 9(n) _ 1 (3) 
C n+1 

It is easy to check that e n ~ 1/n for an asymptotic series of the form (|l|). When the asymptotic 
behaviour of e n is known, it is possible to write down an asymptotic formula for the relative 
error S[n/m\ in the PAP estimate c^f +M+1 of the perturbative coefficient cjv+m+i- 

_ d$+ M+ i — cn+m+i 

0[N/M\ = (4 J 

CN+M+1 

When e n ~ 1/n, as is the case for an asymptotic series of the form ([I]), we have been able to 
demonstrate, for all values of and several values of M, that 

M! 

S [N/M ] -~Tm' where L = N + M + a'M + b . (5) 
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for some choice of numbers a', b. This implies that asymptotically 



hi | 5[ M /m] I — —M [1 + ln(2 + a')] (6) 

We have verified that this formula is numerically accurate for simple series of the form ([!]) 
in which c n = n\K n rC . Moreover, we have checked the prediction (|6]) for several different 
asymptotic series, including that for the QCD vacuum polarization D function in the large- Nf 
approximation || , where it agrees numerically very well with the relative error reported in panel 
(a) of the figure in Ref. We note that the large- Nf D function contains an infinite number 
of renormalon poles, which could in general provide important corrections to the leading-order 
formulae (f|, |]). The fact that this is not the case supports the empirical utility of the PAP's 
even beyond the idealized analytical case ([[]). 

In the previous paragraph, we have discussed the use of Pade approximants to estimate 
the next term in a given perturbative series, and, as we have seen, sufficient conditions for 
the convergence of such PAP's are known, so that the main open issue is their actual rate of 
convergence in practical applications. Next we discuss how to use PA's to estimate the "sum" 
of the perturbation series, which we term Pade Summation (PS). 

It is important to note that convergence of the PAP's is largely independent of the "summa- 
bility" of a given series. Thus, for example, the PAP method gives equally precise predictions 

oo oo 

for the next terms in the series ^n\x n and J^n! (— x) n , even though the latter is Borel 

o o 

summable, while the former is not. 



The formal issues related to Pade Summation are less clear, especially because most per- 
turbation series of practical interest are not Borel summable (since their Borel transforms have 
poles on the positive real axis). One well-defined prescription for defining the "sum" of such 



a series is the Cauchy principal value of the inverse Borel transform integral ||10|| , so we ask 
whether PS "converges" to this prescription for the "sum" of the series. Such series are in 
general obtainable from functions with cuts on the positive real axis, with a toy example being 

oo 

provided by the simple asymptotic series Y^n\ x n , which is a formal expansion of 

o 

e -t 1 ,oo e -v/x 

dt=- aly 7 

1 — xt x Jo 1 — y 

We see from the representation on the right-hand side of (|7]) that this series corresponds in 
QCD language to a single simple infrared renormalon pole, whilst the left-hand side of (0) 
exhibits a cut on the positive real axis. Figure 1 exhibits the errors with respect to the Cauchy 
principal value of the integral (|7]) of conventional partial sums and PA's to the series Sg°n!x n 
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for the case x = 0.1 for which n opt = 9. We see that the first few PA's have relative errors that 
are considerably smaller than those of the partial sums. The best PA has an error comparable 
to An opt , despite using as input a number of terms that is less than n opt . However, we see that 
higher-order PA's exhibit relative errors that are irregular, and may even be less accurate than 
the conventional partial sums. This is because the Pade method mimics the cut on the left- 
hand side of (0) with an ever-denser set of poles, the [N/M] PA's having typical separations 
Ax = 0(1/ M). This makes it ever more difficult to avoid nearby poles when the PA's are 
evaluated at any fixed value of x, and therefore to define "convergence" of the PS procedure. 
However, PA's of sufficiently high order are again much more accurate than the partial sums, 
which are blowing up because of the renormalon singularity. 

As also seen in Figure 1, it is possible to improve on the simple PA's. One way of smoothing 
out the irregularities associated with nearby poles is to evaluate the real parts of the [N/M] 
PA's off the real axis[] at x+ie. It is a theorem |nj that this smoothed PS prescription converges 
to the Cauchy principal value in the cut plane | arg z\ > 0, if the weight in the integral over the 
positive real axis is positive. Another approach is simply to remove from the PA any spurious 
contribution of a nearby pole. Yet another is to make a Taylor expansion of the PA and truncate 
it at order n opt . We have developed |L2| criteria for deciding which of these is better for any 



given case, and the application to our toy example is also shown in Figure 1, designated by 
"combined method" . Since the investigation of these techniques is still in progress, we do not 
discuss them further in this paper, deferring this to a future publication [T2[], where a detailed 
description will be given, together with physical applications. 



3 Application to the Bjorken Sum Rule 

We now discuss a concrete application of PA's to a perturbative QCD series, namely that for 
the Bjorken sum rule, which takes the following form in the MS renormalization prescription 

0,1: 

f\g p l (x 1 Q 2 )-g n l (x,Q 2 )]= l -\g A \f(x) : f(x) = 1 — x — 3.58x 2 — 20.22x 3 + ... + (HT) (8) 
Jo o 

for Nf = 3, as relevant to the Q 2 range of current experiments, where x = a s (Q 2 )/ir, the dots 
represent uncalculated higher orders of perturbation theory, and (HT) denotes higher-twist 
terms. The perturbative series in (|8]) is expected to be dominated by renormalons in large 
orders |TJ]], leading to growth in the perturbative coefficients of the form shown in ([!]). The 
1 We thank A. A. Migdal for pointing this out. 
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PA's to the series @ yield the following predictions for the next term 

Z 4[PA] 



4kl«-lU ([1/2] PA) 



(9) 

([2/1] PA) 



in this series. 



In order for these PAP's to be useful, it is important to estimate the errors involved. The 
asymptotic error estimate given in eq. ([5]) requires as input values of a' and b, which are not 
known a priori. Experience with many series shows that typically and b « 0. 

With a' = and b = we obtain a ballpark estimate 

6 m ^ -2/9; 6 m ~ -1/3 (10) 

In a previous paper we used a different method to estimate the errors of the Pade predic- 
tion @, obtaining cf&> A i = —112 ± 33 as the error- weighted average of the [1/2] and [2/1] 
approximants. We also pointed out that this prediction is close to an estimate made using the 
Effective Charge Method (ECH): 

We now note further that even though the error estimates ( |T0"D are in principle expected to 
hold only asymptotically, in practice (§) and ( |i~l~l) are consistent with each other within these 
estimates. We take this as an indication that the true value of cf J is likely to be in the range 
predicted by the PAP and ECH methods. Moreover, as we shall show in section 4, it seems 
that the perturbation series in (|j) is already dominated by a single infrared renormalon, in 
which case the quantities (|i~0"D are no longer "statistical" errors, but fractional corrections to 



be subtracted from (g), improving the concordance with the ECH estimate (11 



The next step is to apply the PS procedure to estimate the complete correction function 
f(x) in (H). Figure 2 compares /' 2//2 l(x), the [2/2] PS estimate of f(x) (obtained using the ECH 
value QllD of the fourth-order perturbative coefficient), 

, [2/2l , v _ 1 - 8.805 x + 11.974 x 2 

J W- i _ 7. 8 05 x + 7.753 x 2 { ' 

with the [1/2], [2/1] PS's and with the partial sums of the perturbative series up to order x 3 
and x A (the latter also taken from (|TT|) ). We see that the different PS's are numerically quite 
stable in the range x ^ 0.1 of relevance to present experiments, which is related to the fact that 
the nearest poles are some distance away (x = 0.18 for the [2/1] PA, x = —3.41 and x = 0.18 
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for the [1/2] PA, and x = 0.15 and x = 0.86 for the [2/2] PA). This means that the "combined 
method" for smoothing of PA's described at the end of the previous section is not necessary. 

We now check the reliability of the PS estimates of f(x) in two different ways, first by 
checking their renormalization-scale dependences. Figure 3 compares the values of a s (Q 2 ) 
estimated from a fixed value of f(x) at Q 2 = 3 GeV 2 , when one renormalizes at scales \x in 
the range between Q/2 and 2Q, using the PS or the third- or fourth-order partial sum. We 
see that the value of a s (Q 2 ) extracted from the PS is indeed much less /i-dependent than the 



values extracted using the partial sums |15|], consistent with our belief that the PS's provides 
a reliable estimate of the full correction factor f(x), which should be independent of fi. In 
particular, the /i-dependence of the [2/2] PS is small compared with other possible sources of 
theoretical error. 

We have also compared the PS's with results based on the BLM treatment of the perturba- 
tive series in the MS scheme, in which the growing higher-order coefficients are absorbed into 
the scales Q 2 at which x = a s {Q 2 )/^ is evaluated in each of the lower-order terms. As de- 
scribed elsewhere [|12| , we find that our PS procedure agrees very well with the BLM procedure 



when applied to the perturbative series for the Bjorken sum rule, adding further support to our 
contention that the PA's may indeed accelerate usefully the convergence of perturbative QCD 
series, as suggested by the general arguments of Section 2. 



4 Pade Approximants in the Borel Plane 

If one knows the asymptotic behaviour of the series under study, one can go further. In 
particular, if the perturbative coefficients diverge as in (|I]), which is believed to be the case in 
perturbative QCD, corresponding to a discrete set of renormalons 0, it is useful to consider 
the Borel transform of the series S(x) in eq. (Q): 

00 r { A\ n+1 

s( y ) = J2 tnv n ■ £n = hr ; A = (33 - 2 jv»/3 (is) 

n=0 n - V^O/ 

If the Borel transform S(y) indeed has a discrete set of renormalon singularities rk/(y — yk) P , 
where the rjt's are the residues, PA's in the Borel plane are a priori well suited to find them. 
Indeed, if there is a finite set of renormalon singularities, as occurs for the Bjorken sum rule 
series in the large-iV/ approximation, higher-order PA's will be exact. In general, the removal 
of the n\ factors in the coefficients (|T3|) means that the corresponding quantity measuring the 
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rate of convergence of the Pade prediction for the next term is 

C n C n +2 1 



C n+1 n2 



(14) 



which is much smaller than the previous e n ~ 1/n (|3]). This means that the relative error <W,mi 
in the PA of the Borel-transformed series is also much smaller than (pj) for the original series, 



(Ml) 
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5 [M/M ] ^ - -jjm- (15) 

corresponding asymptotically to 

In | 5 [M /M] | ^ -2M[1 + ln(2 + a')} (16) 

We have also checked this prediction for several different asymptotic series, including that 
for the QCD vacuum polarization D function in the large- Nf approximation ||, which has a 
discrete infinity of renormalon poles. The prediction (^) again agrees numerically very well 
with the relative error reported in panel (b) of the figure in Ref . M , which is much smaller than 
that for the naive PA in panel (a). Again, the success of the prediction (|l5|) gains significance 
from the fact that the large- Nf calculation exhibits an infinity of renormalon poles, indicating 
that the Borel PA's are useful in the real world, and not only in idealized simplified situations. 

We now apply this combined Borel/Pade technique to the QCD perturbation series (|8]) for 
the Bjorken sum rule. With the normalization of the Borel variable y implicitly defined through 
eq. flTB|), the [2/1] PA to the Borel transform of (|§|) has a pole at y = 1.05 with residue r = 
0.98, as seen in Fig. 4. The appearance of a pole near y — 1 is encouragingly consistent with 
the exact large-iV/ calculations [13], which yield poles at y = ±1,±2. In the MS prescription 



that we are using, the residues of these poles contain factors exp[5y/3]. Therefore, it is not 
surprising that an infra-red renormalon pole at y = 1 emerges more clearly than an ultraviolet 
renormalon pole at y = —1. Prima facie, the message of this analysis is that the calculated 
Bjorken series is already dominated by the expected leading infra-red renormalon. 

Encouraged by this success, we have made fits to the Borel transform of the Bjorken series 
with varying numbers of poles whose locations are fixed in accordance with theoretical expec- 
tations. The residues found in these various fits are also plotted in Fig. 4, where we see the 
following points: 

(i) the residue of the y = 1 pole is consistently found to be positive and around unity, 

(ii) the residue of the y — — 1 pole is much smaller, and consistent with zero, 
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(iii) there is room for a second pole at y = 2, but it is not possible to disentangle this from a 
higher- lying pole. 



We now discuss the possible implications of this Borel/Pade exercise for phenomenology. 
First, we note that the dominance by a single infrared renormalon pole indicates that, as already 
remarked, the fractional errors ( |I"0|) should be subtracted horn the naive estimates (||), bringing 
them into better agreement with the ECH estimate (|TT|): 

<& PA] --136,-152 (17) 

from the [1/2] and [2/1] PA's respectively. Secondly, it is well known that the magnitude of the 
residue r 1; of the y = 1 pole corresponds to a possible renormalon ambiguity ±7rri relative to 
the Cauchy principal value discussed previously. This is also shown for our toy series in Fig. 1. 
Taking r\ for the Bjorken series from Fig. 4, we find an ambiguity 

A(r?-r?) = ±^0.98vr^ (18) 
6 Q z 

in the perturbative contribution to the Bjorken sum rule. Numerically, with A = 250 ±50 MeV, 
this corresponds to ±0.040 ±0.016 GeV 2 /Q 2 , which is to be compared with previous QCD sum 
rule higher-twist estimates that yield pi 



(see also the discussion pertaining to eq. (7) in Ref. |I7[], and a recent estimate in Ref. |H 
We note that the order of magnitude of the renormalon ambiguity we find (|T8| ) is close to the 
higher-twist calculation flI5|). 

Since the full QCD prediction for any physical quantity must be unique, the renormalon 
ambiguity must be cancelled by a corresponding ambiguity in the definition of the higher-twist 
term. For this reason we do not interpret the renormalon ambiguity as leading directly to an 
ambiguity in a s , but rather use the uncertainty in the higher twist term (|T9|). 



5 Extraction of a s from Bjorken Sum Rule Data 

We conclude this paper by extracting ]17| a s (M§) from data on the Bjorken sum rule, including 
a discussion of theoretical errors. We combine the available experimental evaluations of Y p { n (Q 2 ) 
to obtain 

r?(3GeV 2 ) - r™(3GeV 2 ) = 0.164 ± 0.011 (20) 
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where we have evolved all the experimental results [19]- [21] to the common reference scale 



Q = 3 GeV . Evolution of the proton data from 10 to 3 GeV is based on the analysis 



in Ref. p0|| . The corresponding evolution of the deuteron data is based on the proton data, 
combined with the expected Q 2 dependence of the Bjorken sum rule, with a s (Q 2 ) needed as 
input determined via iteration, in a self-consistent way. Evolution to 3 GeV 2 from values of 
Q 2 different from 10 GeV 2 was treated in a linear approximation to the full dependence of the 
data on 1/ log(Q 2 ), which is adequate for this purpose. 

The range of f(x) corresponding to ( |2C|) is shown as a vertical error bar in Figure 2. We use 
the [2/2] PS (HD to obtain: a s (3 GeV 2 ) = 0.328±g;g|f. We evolve this up to M| by numerical 
integration of the three- loop (3 function ||, locating the 6-quark threshold in the MS scheme 
at nib = 4.3 ± 0.2 GeV ]^2][ , and using the three-loop matching conditions of ref. |23|, to find 



a s {M 2 ) =0.1191°;^ ± ••■ , (21) 
where the ± . . . in (EI]) recalls that theoretical errors remain to be assigned. 



There is a theoretical error associated with the spread in the different evaluation procedures 
shown in Fig. 2, which we estimate from the difference between the [2/2] and [1/2], [2/1] PS's 
to be A proc a; s (3 GeV 2 ) = ±0.014, corresponding to 

A proc a s (M 2 z ) = ±0.002 (22) 

Another way of looking at the theoretical error uses the /i-dependence of the [2/2] PS shown 
in Fig. 3 to estimate A^a s (3 GeV 2 ) = ±0.009, corresponding to 

A M a s (Mf ) = ±0.001 (23) 

from varying [i between Q/2 and 2Q. Both ( p2|) and (|23|) are estimates of the uncertainty in 
(Mf ) due to our uncertainty in the functional form of the QCD correction factor f(x), so that 
it could be regarded as double counting to include them both. Nevertheless, to be conservative 
we will add them in quadrature.^] 

We also include a shift and error in the determination of a s (Q 2 ) inferred from the estimated 
range ( |T9"D of the higher-twist correction: A#t«s(3 GeV 2 ) = —0.024 ± 0.014, corresponding to 

A HT a s (Ml) = -0.003 ± 0.002 (24) 



2 We have also considered possible systematic theoretical errors due to uncertainties in the evolution of a s 
up to Mz-, including unknown higher-order terms in the QCD /3-function, the uncertainty in mj, and freedom 
in treating the heavy-flavour threshold |p3[ . They contribute an error in ct s (Mz,2) which is much less than 
±0.001 . 
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Combining 



we extract 



a 



0.116±g:gg| ±0.003 



(25) 



where the first errors are the experimental errors in 



and the second errors are the sums 



in quadrature of the errors in (22|), (|23), (j24|) . 



Our final value of a s (M|) (|25"1) is compatible with the central value extracted from compi- 
lations of previous measurements |24J], and has an error which is competitive. As well as the 
experimental error, we have included motivated estimates of a number of theoretical errors, 
using information obtained from our study of Pade approximants. We believe that this exercise 
demonstrates the utility of Pade Approximants in QCD, and the value of polarized structure 
function data for determining a s (M|). 
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FIGURE CAPTIONS 



Fig. 1. The relative errors between partial sums of the series S(x) = T,n\x n and the 
Cauchy principal value of the series (solid line) is compared with the relative errors of Pade 
Sums (dotted line). We see that the relative errors of the Pade Sums are smaller than those of 
the partial sums in low orders, fluctuate in an intermediate regime, and are again more accurate 
than the partial sums in higher orders. The fluctuations are associated with nearby poles in 
the Pade Sums, that may be treated by the "combined method" mentioned in the text, shown 
as the dashed line. 

Fig. 2. Different approximations to the Bjorken sum rule correction factor f(x), third- 
order and fourth-order perturbation theory, [1/2], [2/1] and [2/2] Pade Sums ([12]), are compared. 
Also shown as a vertical error bar is the value of f{x) we extract from the available polarized 
structure data ([201). 

Fig. 3. The scale dependence of a s (3GeV 2 ) obtained from a fixed value fix) = (G/gA) x 
0.164 = 0.783, (cf. eq. (|20|) ), for Q/2 < \i < 2Q, using the naive third- and fourth-order 
perturbative series and the [1/2], [2/1] and [2/2] PS's. 

Fig. 4. The locations and residues of poles in the [2/1] PA and in rational-function fits to 
the Borel transform of the first four terms in the perturbation series for the Bjorken sum rule. 
We note that the location of the lowest-lying infrared renormalon pole is estimated accurately 
by Pade Approximants in the Borel plane, and that its residue is stable in the different fits. 
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